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Abstract 

In this paper we present a formulation of the nonlinear stochastic differen- 
tial equation which allows for systematic approximations. The method is not 
restricted to the asymptotic, i.e., stationary, regime but can be applied to de- 
rive effective equations describing the relaxation of the system from arbitrary 
initial conditions. The basic idea is to reduce the nonlinear Langevin equa- 
tion to an equivalent equilibrium problem, which can then be studied with 
the methods of conventional equilibrium statistical field theory. A particular 
well suited perturbative scheme is that developed in quantum field theory by 
Cornwall, Jackiw and Tomboulis. We apply this method to the study of N 
component Ginzburg-Landau equation in zero spatial dimension. In the limit 
of N — > oo we can solve the effective equations and obtain closed forms for the 
time evolutions of the average field and of the two-time connected correlation 



functions. 
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I. INTRODUCTION 



Non linear stochastic differential equations (SDE) in one or more variables occur fre- 
quently in the description of a wide variety of physical phenomena such as those involving 
relaxation towards a steady state |I],0]. A celebrated example of SDE is the Langevin equa- 
tion, where one relates the rate of change of some physical observables to a drift term, i.e. to 
a deterministic driving force, plus a stochastic noise. Whereas stochastic linear differential 
equations are in principle amenable to analytic solutions, those of nonlinear type are much 
more difficult to treat and in most cases one has to resort either to computer simulations 
or to approximate schemes. Several approximations exist in the literature for the analysis 
of nonlinear SDE. A wide used method is the so called statistical linearization [§] which 
consists of the replacement of the nonlinear SDE by equivalent linear ones whose coefficients 
are determined by some error minimization algorithm. This method in its original form is 
appropriate in the asymptotic, i.e., stationary, regime. In many physical systems, however, 
it is the transient regime in which one is interested. To deal with this problem the idea 
of statistical linearization has been extended to include a possible time dependence in the 
parameters . The drawback of these methods is that they are not derived in a systematic 
way so it is not simple to improve them. 

An alternative procedure is the so called dynamical Hartree approximation. This scheme 
amounts to neglect correlations of higher order than the second and is exact whenever the 
probability distribution associated with the problem is Gaussian. This is not the case of 
strongly interacting systems, so corrections to the Hartree approximation should be taken 
into account. Also in this case it is not simple to improve the approximation. 

In this paper we present a formulation of the nonlinear SDE which allows for system- 
atic approximations. This is achieved by reducing the nonlinear Langevin equation to an 
equivalent equilibrium problem, which can be analyzed with the methods of conventional 
equilibrium statistical field theory. In particular we have applied a method originally de- 
veloped in quantum field theory by Cornwall, Jackiw and Toumboulis ||[7|], alternative to 
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conventional perturbation theory, because a normal coupling constant expansion can only 
be used for the study of small corrections to the deterministic result. In this respect the 
present approach is alternative to the field theoretical treatment based on the introduction 
of fermionic auxiliary fields. 

To illustrate the method we shall study an N component Ginzburg - Landau equation 
in zero spatial dimension with the purpose of deriving systematically the time dependent 
Hartree equations and the first corrections. The same model has been discussed previously 
by Scalapino and coworkers |§, who introduced an approximation scheme for the Langevin 
equation based on the expansion in the small parameter 1/N. We shall obtain a solution 
which represents a systematic treatment, and as a bonus we derive a closed equation for 
the two-time connected correlation function. This will be discussed in Sects. |T|, [TV] and [V], 
where the corrections to the N — > oo solution are considered. 

The formalism is established in Sec. |I| where we construct the generating functional 
T for the average value of the observables and its correlations. In Sect. |VT| we provide 
an alternative derivation of the Hartree equations, based on a variational principle similar 
to the Feynman method in equilibrium statistical mechanics, and discuss the relevant case 
N =1. 

II. THE FORMALISM 

In this Section we derive the effective action formalism for the Langevin equation. The 
path integral method constitutes a convenient representation of the Langevin equation for 
a field 4>(t). Within this approach, the original stochastic differential equation, where the 
depends on another field £, called the noise, is reformulated by constructing an effective 
action for the field only obtained by eliminating the noise. The advantage of this trans- 
formation is that one can employ the well known methods of equilibrium statistical field 
theory. To keep the notation as simple as possible the derivation will be carried out for a 
single component real field. The extension to iV-vector fields will be discussed later. 
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The time evolution of the field </>(£) is governed by the Langevin equation: 

9 ^ 9 QUI . C 

where S[<p} is an "energy" function, and £ a Gaussian random variable with 



(£(*)> = o, <e(*)e(0> = r*(t-f)- 



(i) 



(2) 



We shall derive a generating functional from which the correlations can be obtained. Pro- 
ceeding in the standard way, see e.g. Ref. 0, we introduce an external source J and define 
the generating functional 



z[j] = N J v"<pvsv[<KQ)W-<f>e) 

- r dt 

. Jo 



x exp 



cxp 



e 2 

dt — 

2T 



(3) 



where 0£ is the solution of stochastic eq. ([TJ) subject to some set of initial value conditions 
0(0) assigned with probability P[0(O)], and J\f is a normalizing constant. The functional 
integral on in eq. (f|) includes integration over 0(0) and 0(r). We denote it by the double 
quote: T>"<p. The 5-function stands for 



5(0-0 5 ) = 5 



d<p dS _ 

— H f 

dt 86 ^ 



det 



(4) 



where det |<5£/<50| is the Jacobian of the transformation £ — > 0. With well-known manipula- 
tions, see e.g. Ref. one nas 



det 



cxp 



1 r a 92 s 

2 Jo w 



(5) 



In deriving eq. @ we have used the forward time propagation Green function 9(t — t') of 
the operator dt. Moreover we have used the definition 



m = 5 . 

This choice corresponds to the "physical" regularization of the noise term £ as 



(6) 



(?) 



where rj(t) is an even function sharply peaked at t = 0, whose integral from — oo to +00 is 
equal to 1. The 5-correlated noise is obtained in the limit of vanishing width. In terms of 
stochastic differential equations this corresponds to the Stratonovich formalism [ |TTf . 

At this stage one eliminates the noise field by inserting eq. (Q) into eq. (Q) and performing 
the £ integral over the noise obtaining, 



Z[J\ = N / £>"0P[0(O)] exp 



1/90 dS_\ 

> zr{dt + dj) 

1 r , d 2 6 r , • 



(8) 



The argument of the exponential can be simplified by performing the integration of the term 
Jq dt<j)dS/d(f) = S[(/)(t)] - 3[<f>(0)], so we finally have 

Z[J]=Af f P0(O)P[0(O)]e 5 ^/ 2r P0(r)e 5 ^/ 2r 



x X>0 exp 



-IU) - dtJ<f) 
Jo 



(9) 



where the action I(<ft) is given by 



.dt + [d6, 



1 r 

2 7o 90 2 ' 



(10) 



and denotes integration over all paths starting at 0(0) for t = and ending at 0(r) for 
t = r. It is defined as 

N-l 

N^oo 



V(p= lim Yl d</>(ti 



(11) 



i=l 



where <f>(ti) is the field at time t% = ie, having sliced the interval to r in iV parts of size 
e = t/N. 

In eq.(^) the integration over the end points just fixes the boundary conditions at t = 
and t = t, so without lost of generality, we can consider a 'reduced' generating functional 



Z[J] = M T' 01 £>0 exp 



dt J0 



(12) 



where now the integral runs over all paths which start at t = from 0(0) = 0o and end 
at t = r at 0(r) = 0i. In the limit r — > 00, or equivalently t/r 1, the path becomes 
independent of the final value 0(r). 
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We note however that the presence of the additional constraint is necessary to select 
paths which are solution of the original equation of motion eq. ([!]). Equation ([H]) and flT2|) 
lead in fact to second order differential equations of motion, whereas the original stochastic 
equation is of the first order. The additional constraint at t — r makes the problem well 
defined since the paths in eq. ( [T2| ) must satisfy the two constraints 0(0) = 0o an d 0(t) = <f>\. 
Once the two boundaries conditions are imposed the path is also solution of the first order 
differential equation ([]]), as can be easily seen in the limit r — > 0, i.e., the deterministic 
limit. 



In general the calculation of path integrals like eqs. (|T0|) and (]12|) is not at all straight- 
forward. Nevertheless quantities of physical interest can be obtained. In our case we are 
interested into the noise-averaged value of the field (4>{t)) and correlations (<f){t) <p{t r )) as a 
functions of time. An advantage of the present formalism is that self-consistent, systematic 
variational principles for these quantities can be obtained using a method introduced by 
Cornwall, Jackiw and Tomboulis in Quantum Field Theory ||. The basic idea is to derive 
an effective action which is stationary at the physical values of {4>{t)) and {4>(t) 4>{t')). 

The method starts by generalizing eq. flE2| ) to account for the composite operator 
<p(t)(f)(t'). We then define the generating functional 

Z[J, K]=M V(f) exp - dt J(t) (pit) 

Jo,^ 1 JO 

- \ r dt / T dt'0(t)ir(t,t')0(o] (is) 

Z Jo Jo J 

where J and K local and a bilocal source, respectively. 

By taking functional derivatives with respect to the external sources the averaged corre- 
lations of (p can be obtained. In particular, by considering W[J, K] — — In Z[J, K], we have 
for < s < t 

jj^WllK] =<0(5)>=g(5) 

(14) 

j^-^ W[J, K] = i(0( S ) 0( S ')> = | [q(s) q(s') + G(s, s')] 



where the averages are obtained with the weight of eq. (]T3|). In the limit of vanishing 
external sources q and G become the noise-averaged field (<f)(s)) and connected two-point 
correlation function (<f>(s) <f>(s r )) c of the process described by the Langevin equation (|I|). 
By Legendre transforming W[J, K] we can eliminate J and K in favor of q and G: 

r[q,G] = W[J,K}- f dsq(s)J(s)-l f ds f ds' q(s) K(s, s') q(s') 
Jo 2 Jo Jo 

-- [ T ds [ T ds'G{s,s')K{s,s') (15) 
2 Jo Jo 

where J and K are eliminated as a function of q and G by the use of eq. ( fH] ) . It can be shown 
that F[q, G] is the generating function of 2PI (two-particle irreducible) Green functions, i.e. 
it is given by all diagrams that cannot be separated in two pieces by cutting two lines 
The external sources can be obtained from T[q, G] as 

() ' T[q,G] = -J(s)-f T dsK(s,s')q(s) 



Sq(s) 

(16) 
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r[g,G] =-\K(a,s!). 



[ 5G(s, s ) J - 2 



The physical process corresponds to vanishing sources, J = K = 0. From eq. (0) it 
follows that in this limit the value of q and G are determined by the stationary point of 
F[q, G}. We have thus obtained a variational principle for the the noise-averaged field (0(s)) 
and connected two-point correlation function (0(s) <f>(s')) c of the process described by the 
Langevin equation QTJ) . 

The next step is to evaluate T[q, G}. Following Ref. PJ^] T[q, G] can be written as 

T[q, G] = I(q) + i Tr In G' 1 + i Tr P" 1 ^) G + r 2 (g, G) + const. (17) 



where 7(g) is given by eq. (|T0|) with — > g, 



<J0(s)50(s') 



(18) 



50( S )50( S ') 

with propagator of the "free" theory. The functional T 2 is given by the sum of all 2PI 
vacuum diagrams of a theory with interactions determined by L mt and propagators G. The 
interaction term is defined by the shifted action 
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I(q + <P)-I(q)-<P 




(19) 



5<f) 



<p=q 



This procedure corresponds to a dressed loop expansion with vertices which depend on 0, 
and can thus exhibit non-perturbative effects even for a small number of dressed loops. 
The crucial point is that it in no sense corresponds to a perturbation theory in physical 
amplitudes. The physics is obtained by going to the stationary point of the expansion with 
respect to both variables <fi and G. This leads to nonlinear dynamical equations for and 
G. If one could sum up the whole series, the exact value of <fi and G will emerge from the 
stationary point. If the series is truncated one gets approximate values of <fi and G. The 
bonus is that on can get equations describing nonperturbative behaviors. The drawback is 
that in general there is no systematic knowledge about errors occurred. 

In the next section we apply the above formalism to a O(N) problem where the leading 
contribution to I^ can be extracted in the N — > oo limit. 



To illustrate the formalism introduced in the previous section we consider an N- 
component Ginzburg-Landau time-dependent field fa with quadratic local interaction in 
zero spatial dimension. When discussing fluctuations effects to any given order in a pertur- 
bation expansion one is not usually able to justify the neglect of yet higher orders. However 
for theories with large N internal symmetry group there exists another perturbative scheme, 
the 1/N expansion. The model is specified by the evolution equation 



III. THE MODEL 



d ^ 9 Q\M -L- P 

-m=~dfa m + 



(20) 



(21) 



where Gaussian random noise is defined by: 



<&(*)> = o, 



(22) 
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and we assume a < and A > 0. Generalizing eq. (|H]) to iV-component fields one finds the 
following action 



J(0) = i f T dt 

1 Jo 



4> , m2 ^2, X (rk 2\2 i 90 ( ,2s 3 aN 



where the parameters m 2 , Aq and go are related to the original constants by: 



2 2 

m = a — 



AT AT 

6 ~ 3iV' 



(23) 



(24a) 



A = 4 aA, 



(24b) 



p = 10 A 2 . 



(24c) 



The last term in eq. ( p3]) does not depend on and can be absorbed into the definition of 
the normalizing constant J\f in eq. (f0f). 

In the limit iV — > oo we can calculate explicitly the leading order term of the functional 
(|l~7|) following the same steps of Ref. ||12|1 . The "action" ( p3l) corresponds to a classical 
6 -theory in one spatial dimension. 

From eq. (p3|) it follows that the leading contributions to eq. ( |17|) for iV — > oo are 



? 4 (*) 



G(t,t')5(t-t') (25) 



and 



m + — g(t) + n-^g (t) 



5(t-t')- 



dt 2 3!i\^ v ' 5!iV 2 

The leading order 2PI diagrams iV — > oo, shown in Fig. [3], lead to 

r 2 (g,G) = ^ f diG 2 (M) + ^f fdtG 3 (M) + ^ [ T dtq\t)G 2 (t,t) 



(26) 



4! Jo 



6! Jo 



6! Jo 



(27) 



where q(t) = (<pi(t)) assuming that the symmetry is broken along the direction '1'. 

Stationarity of the functional T[q, G] with respect to q(t) and G(t, t') yields the dynamical 
equations for the order parameter and its fluctuations which read, respectively 
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d 2 



dt* +m+ 3!iV 



^0 2/jA i 90 4/,\ 



5!iV 2 



+ ^(M) + |^(M) + 5 , v 



2 % 2 WG(M) 



?(*) = o 



(28) 



and 



~W +m + 3\N 
Ar 



5!iV 2 



+ i t} + i G2(t; t} + If g2(t) t} 



G(t,t') = r<y(t-t') 



(29) 



These coupled dynamical equations are exact to leading order in N. 

The effective dynamical equations ( p8|) and (p9|) can be also derived by a variational 
approach to the path integral (||), or Hartree approximation, where one seeks for the best 
quadratic approximation for the action. 



IV. THE SOLUTION 



To solve the coupled effective dynamical equations (p8|) and (|29|) we rewrite them as 
d 2 q(t) 



dt 2 
d 2 G(t,t') 

dt 2 

where 



F[q(t),G(t,t)]q(t) (30) 
F[q(t), G(t, t)] G{t, if) - T8{t - f) (31) 



F[q(t),G(t, t)} = m * + A ,2 (t) + ^) 

+ ^ G(t > t] + 1 G2(t> t] + H g2(t) t} (32) 

We are eventually interested into solutions for times t such that t/r — ► 0. Under this 
assumption the effective dynamical equations (|30"D and fl3~I|) can be reduced to simpler first 

order non linear differential equations by using the following representation for q(t) and 
G{t,f): 

q{t)=q(0)fi(t) (33) 
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and 



G{t, f) = h(f) f 2 (t) 9(t' -t) + h(t) f 2 (t') 6{t - t') (34) 

with 

/ 1 (t)=e-/o^W (35) 

/ 2 (t) = re-/oe 2 i T ^ (36) 



where the function R(t) is solution of the first order non linear differential equation 

R 2 (t)-F[q(t),G(t,t)] (37) 



Of 



By inspection of eqs. (|32| ) and ( p7|) it follows that R{t) should have the functional form 

R(t) = aC(t) +{3q 2 (t) +7 (38) 

where C(t) = \im t ^t> G(t, t'). The parameters a, [3 and 7 are determined by substituting 
R(t) from eq. (|3"8"| ) into eq. (P^), and eliminating dC(t)/dt and dg(t)/dt with the help of 

MS = -J2(t)g(t) (39) 
^ = -2i?(t)C(t) + r (40) 

obtained from eqs. ( p3| ) - (|36|) . We obtain the following solution for R(t) 

R(t) = ±C(t) + ± q 2 (t) + a. (41) 

In principle there exists another set of parameters, but it leads to an asymptotically unsteady 
state. 

The first order non linear differential equations (|39D - (f|l]) gives the full description of 



the model (^) - fl2~l|) in the limit iV — > 00 for all times t such that t/r — > 0. 

If q(t) is not identically equal to zero it is not straightforward to solve analytically the 
set of equations (^) - fl41|). Nevertheless these can be easily solved numerically for any set 
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of initial conditions. We note that this is not the case for eqs. ( pq) and (|29|). These, indeed, 
suffer of strong numerical instability and one has to resort clever algorithm to handle them. 

On the contrary if q(t) = for all times, then we can find a closed analytical solution. We 
note that from the structure of the equations on can see that the O(N) symmetry dictates 
the equilibrium value of q{t): lim t ^ 00 g(t) = 0. Thus if we assume that q(0) = then the 
solution is q(t) = 0. In this case the equation for G(t, t') can also be solved, and if we take 
as initial condition C(t = 0) = 0, the solution reads 



fi{t) = e- at/2 

3r 



1 - (a//3) 



\ e XAt/3 _ ( a /0) e -\At/3 



(42) 



at/2 



A A \ e AA '/ 3 - (a//3)e- AAi / 3 



n/ h sinh(AAi/3) (43) 



where 



3a / 3T 3a / 3T 

a = ~T + V + T = ---t' + - (44) 



Substitution of eqs. (|4|) - (gj) into eq. (|34[) leads to the solution for G(t,t'). In the limit 
£' — > t we recover the result of Ref . || for C (t) . 

V. BEYOND THE HARTREE APPROXIMATION 

The next leading terms of order 1/N can be included systematically by evaluating dia- 
grams not included in eq. ([27|). In the case q(t) =0 the work is simplified since one does not 
have to consider separately transverse and parallel components of the correlation function 

(0iOW))c 

The diagrams contributing to the first corrections to I^ are shown in Fig. [I] (a) and (b) 
and Fig. ^] and yields 

2 , r „o,. s 6 



T 2 (G) = -A Q J dt G 2 (t, t) + -g J dt G 3 (t, t) 

~\% I dt I dt '° 4 ^ ~ T^^o Jdt J dt'G\t,t')G(t,t)G(t',t' 



Ao gl 
720 



J dt J dt'G\t,t')G{t,t) (45) 
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To improve systematically this result, and the Hartree approximation, one has to consider 
an infinite series of diagrams. A complete summation of the series, see Figs. || and f|, can 
be performed, however, only in the case of a system at equilibrium where time translational 
invariance holds. The final result for T 2 is valid to all orders in A and to first order in 1/N 
and reads O.O]: 



r 2 (G)/r 



iV + 2 



An 



duoi f da; 



4! "i 2n J 2tt 

N + 6 r dull f du>2 f dco>3 
H ^— 9o 



6! 

1 r du n 
+ - / — In 

2 J 2vr 



2tt 



2tt 



2n 



1 + (l + i; G(0 ' 0) ) flM 



where II (a;) is the so called vacuum polarization propagator: 

dtu 



(46) 



2tt 



G(v)G(rj + u) 



(47) 



Upon differentiating with respect to G(u) the functional T[G] with the l/TV corrections 



included one obtains 12 



G- V) = c 2 + m 2 + ^ G(0, 0) + fj G 2 (0, 0) 



+ 



3! 

dcui 
"2tT' 



i + (t + S G (°.°)) n W 



A , go 



G(0,0) -wi) +n(w0 



(48) 



3 30 

Equation ( |48"D represents the spectrum of the equilibrium fluctuations correct to order 1/N. 
The study of these corrections will be the subject of a future publication. 



VI. VARIATIONAL APPROACH 



A well-known approximation often employed in statistical mechanics is obtained by ap- 
plying the Peierls-Feynman-Bogolubov inequality: 



InZ < -lnZ + (C-£o) = W {1) 



(49) 



14 



where £ is action density given by 

I(<P) = I dtC (50) 



o 



and Cq is an arbitrary action density. The average (• ■ -)o in (f|^) is done with respect to 
the probability distribution corresponding to £ an d Z is the partition function associated 
with C . The Hartree method consists of choosing the arbitrary action C to be the most 
general quadratic form: 

f dt£ = 7^ Tdt f T dt' O0i(O - f 'dt (51) 

jo 21 jo jo 77 jo j 

where positive definite kernel. Recalling that: 

W {0) = -\nZ 

= -~ InHti-^OII " £ / T dt f dt'v{t)u-\t,t')v{t') (52) 
2 2 Jo Jo 

Taking the derivatives of with respect to v(t) we find: 

-q(t) = -r / dtV^MXt') (53) 
Jo 



; =-G(t,i') = -r«- 1 (M') (54) 



<fo(t) 5v(t') 
which can be inverted to give: 

[ T dt'u(t,t')q(t') = Tv(t) (55) 

JO 

f T dt' u(t, t') G(t', t") = T5(t- t") (56) 

JO 

The arbitrary functions Uij and V{ are determined by looking for the minimum of to 
have the best estimate of In Z. 

For the case under study by applying the above method to ([H]) one finds 

+ ¥ G(( ' t) + % G ' (i ' i] + fit " 2(i) G(i ' i)] i{t ~ ''>' (57) 

and 
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v{t) = (58) 
Thus inserting (|57|) and fl58|) in eqs. fl55|) and (^) we find the same result as eqs. (|28|) and 



Before concluding the discussion of the Hartree method we shall analyze what happens in 
the case N = 1. Its variational nature in fact justifies its application even when the problem 
under scrutiny does not contain a natural small parameter around which to perform some 
sort of expansion. 

Let us consider the case of a Langevin equation for a scalar field with cubic nonlinearity 
The equations of motion for the average value of the field and for the fluctuations are ||15|| : 



Qf 2 +m (t) + - g (t) 



+ y G(t, t) + f G 2 (t, t) + H q\t) G(t, t)) q(t) = (59) 

[-^2"+ m +Y q + 4! q (t) 

+ y G(t, t) + I G 2 (t, t) + | q\t) G(t, t)} G(t, H) =TS{t- t'). (60) 

With the substitution: 

q(t) = e -fo dTR «V (61) 

and G(t,t') given by eq. (0) we obtain two equations for R q {t) and R{t) analogous to eq. 
(|37l). Using the trial solution 

R q (t) = a q C(t)+P q q 2 (t)+ lq (62) 

and R(t) given by eq.(|38"D we find, up to terms quadratic in the coupling constant A, a q = 
a = A/2, P q = A/6, P = A/2 and •y q = 7 = a. The value of the coefficients coincides with 



the value obtained in the so called Langer-Bar on-Miller approximation |]I6[ to the Langevin 
equation, a result which was also rediscovered few years ago 0], on the basis of a somehow 
ad hoc variational principle. We believe that the present derivation, being based on a path 
integral formulation of the stochastic equations makes the underlying physical assumptions 
more clear. 
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VII. CONCLUSIONS 



Most problems arising in the study of physical problems are most naturally represented in 
terms of systems of nonlinear stochastic differential equations. Many approximation schemes 
have been developed to treat the nonlinear aspect of equations. Usually they are based on 
some reasonable assumptions. The advantage of these approaches is that they may lead to 
relative simple equations. The drawback is that is is not simple to improve the quality of the 
approximation. In this paper we presented an alternative approach to the study of nonlinear 
Langevin equation which allows for systematic development of approximation scheme. The 
basic idea is to reduce the nonlinear Langevin equation to an equivalent equilibrium problem 
to which the methods of conventional field theory can be applied. A particular well suited 
perturbative scheme is that developed in quantum field theory by Cornwall, Jackiw and 
Tomboulis 0. The major advantage is that it leads to a variational principle for the physical 
quantities of interest. 

The method is applied to an N component Ginzburg-Landau equation. In the limit of 
iV — > oo we are able to derive closed forms for the order parameter q(t) and for the two-time 
connected correlation function G(t,t'). We also discuss the first 1/N corrections. The study 
of these is, however, more involved and has not been included in this paper. This is part of 
future work. It will be also of interest to extend the present approach to higher dimensions 
and to explore numerically the predictions of the present approach to finite values of N. 
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FIGURES 
FIG. 1. Leading order 2PI diagrams 

FIG. 2. 2PI diagrams contribution to the first 1/N corrections 

FIG. 3. Three vertices 2PI diagrams contributing to the first 1/N corrections 

FIG. 4. Four vertices 2PI diagrams contributing to the first 1/N corrections 



19 



